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Abstract 

We study the n + 1-dimensional generalization of the dispersionless 
Kadomtsev-Petviashvili (dKP) equation, a universal equation describ- 
ing the propagation of weakly nonlinear, quasi one dimensional waves 
in n+1 dimensions, and arising in several physical contexts, like acous- 
tics, plasma physics and hydrodynamics. For n = 2, this equation is 
integrable, and it has been recently shown to be a prototype model 
equation in the description of the two dimensional wave breaking of 
localized initial data. We construct an exact solution of the n + 1 
dimensional model containing an arbitrary function of one variable, 
corresponding to its parabolic invariance, describing waves, constant 
on their paraboloidal wave front, breaking simultaneously in all points 
of it. Then we use such solution to build a uniform approximation of 
the solution of the Cauchy problem, for small and localized initial 
data, showing that such a small and localized initial data evolving 
according to the n + 1-dimensional dKP equation break, in the long 
time regime, if and only if 1 < n < 3; i.e., in physical space. Such a 
wave breaking takes place, generically, in a point of the paraboloidal 
wave front, and the analytic aspects of it are given explicitly in terms 
of the small initial data. 
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1 Introduction 



The n + 1 dimensional generalization of the dispersionless Kadomtsev - Petvi- 
ashvili equation: 

{ut + uu.^)^ + A±u = 0, u = u{x,y,t), y = (yi, . . . ,yn-i) 



n-1 (^\ 



hereafter referred to as the dKPn equation, describes the propagation of 
weakly nonlinear quasi one dimensional waves in + 1 dimensions. This 
equation was first obtained by Timman [1], for n = 2, in the study of unsteady 
motion in transonic flows, and then derived by Khoklov - Zobolotskaya, in 
the 3+1 dimensional acoustic context in [2], where the main properties of 
the model were pointed out. Equation ([T]) is the x-dispersionless limit of 
another distinguished model, the n + 1 dimensional generalization of the 
Kadomtsev- Petviashvili equation [3], integrable for = 2 by the Inverse 
Spectral Transform (1ST) [U |5] . 

The universal character of ([1]) can be explained as follows. Take any 
system of nonlinear PDEs i) characterized, for example, by nonlinearities of 
hydrodynamic type and ii) whose linear limit, at least in some approximation, 
is described by the wave equation. Then, iii) looking at the propagation of 
quasi - one dimensional waves and iv) neglecting dispersion and dissipation, 
one obtains, at the second order in the proper multiscale expansion, the 
dKPn equation ([T]). Therefore ([T]) arises in several physical contexts, like 
acoustics, plasma physics and hydrodynamics. 

Indeed, consider a nonlinear system of PDEs whose linear limit is the wave 
equation fft' = A/ in n + 1 dimensions; then a localized initial condition 
evolving according to it is concentrated, asymptotically, on the spherical wave 
front: 



x\ -t' = 0(1], \x\ 



\ 



(2) 



k=l 



Looking for nonlinear corrections and quasi - one dimensional propagations 
(say, in the Xi direction), we introduce the convenient variables: 

X = xi — t' = 0(1), yj = e'^Xj+i, j = 1, . . . , n — 1, t = e^'^t', g > 0, (3) 

where e is the order of magnitude of the small amplitude of the wave. Then 
the spherical wave front becomes, approximately, its second order contact. 
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an ellipsoidal paraboloid: 



1 



k=l 



In addition, starting from the dispersion relation of the wave equation and 
looking for quasi one dimensional waves in the Xi direction, we have longer 
wave lengths (smaller wave numbers) in the transversal directions: k± = 
e'^K,±, obtaining that the dispersion relation oj{k) of the wave equation reduces 
to that of the linearized dKPn (II]) (up to the trivial rescaling y — )■ y/\/2): 



y j=i i=2 (5) 

9{x,t') = k ■ X — uj{k)t' ~ kix + K,± ■ y — ji^t. 

This is what happens in the acoustic problem in 3 + 1 dimensions [2], where 
the gas density, the pressure and the component of the velocity in the main 
direction xi are proportional to eu{x,y,t), and u solves the dKP^ equation. 
Starting from the plasma physics equations 

Pf + V-(pw)=0, Vt + {v-V)v + V(f) = 6, A(f)-e^ + p = 0, (6) 

describing a gas of hot electrons on a background of cold ions, where p and 
V are the ion density and velocity and is the electric potential [6], the 
linearized theory leads instead to a generalized fourth order wave equation 
in (3 + 1) dimensions 

A0tt = (Pu - A0, (7) 

reducing to the wave equation in the long wave approximation (small wave 
numbers) : 



1.2 1.2 

^=^lYT¥-k{l-^), e = kl + kl + kl (8) 

Looking, in addition, for quasi one dimensional waves (the wave numbers 
in the transversal directions are smaller than the (small) wave number in the 
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Xi direction: ki = e^Ki, k± = e^~^'^K,±), we obtain 




e{x,t') = k-x-uj(k)t' ~ 



(9) 




where now 



X 



(10) 



If p > 2g, the term /c^ is neghgeable wrt and one obtains again 

the dispersion relation of the hnearized dKP^ equation; if p = 2g, they 
are comparable and one obtains the dispersion relation of the linearized KP 
equation in 3+1 dimensions. 

In the water wave theory, the situation is very similar to that of plasma 
physics and the dKP2 equation is derived in the long wave approximation. 
In all the above three physical contexts, the choice of the exponent q in 
(!3l). (fT0|) comes from the balance, at the second order in the proper multiscale 
expansion, with the nonlinearity of the physical system, and is g = 1/2. 

We remark that the 1 + 1 dimensional version of ([T]) is the celebrated 
Riemann-Hopf equation ut + uux = 0, the prototype model in the description 
of the gradient catastrophe (or wave breaking) of one dimensional waves |7j . 
Therefore a natural question arises: do solutions of dKPn break and, if so, 
is it possible to give an analytic description of such a multidimensional wave 
breaking? 

A first and positive answer in this direction was recently given in fS], for 
the integrable P [IDl IHl 1121 1131 [S] dKP2 case. Indeed, using a novel 1ST for 
vector fields, we have been able to solve the Cauchy problem of dKP2 fT5] 
and of other distinguished integrable Partial Differential Equations (PDEs) 
arising as commutation of multidimensional vector fields [16l [T71 [18]. The 
associated nonlinear Riemann-Hilbert (RH) inverse problem turns out to be 
an efficient tool to study several properties of the solution space of dKP2, like, 
for instance, i) the construction of the longtime behaviour of the solutions of 
the Cauchy problem; ii) the possibility of establishing that localized initial 
profiles evolving according to dKP2 break at finite time and, if small, they 
break in the longtime regime, investigating in an explicit way the analytic 
aspects of such a longtime wave breaking of two-dimensional waves [S]. 
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In this paper, motivated by the analytic results of the integrable case 
n = 2 [8j, we answer the above question in arbitrary dimensions, when the 
model ([1]) is not integrable, under the assumption of small and localized initial 
data. To obtain this result, we first construct an exact solution of equation ([1]) 
containing an arbitrary function of one variable, consequence of the parabolic 
invariance of equation ([T]), describing a wave, constant on its paraboloidal 
wave front, breaking simultaneously in all points of it. Then we use such 
solution to build a uniform approximation of the solution of the Cauchy 
problem, showing that "small and localized initial data evolving according 
to the dKPn equation break, in the long time regime, if and only if 1 < n < 3; 
i.e., in physical space". Such a wave breaking takes place, generically, in a 
point of the paraboloidal wave front, and the analytic aspects of it are given 
explicitly in terms of the initial data. In addition we show that, if the initial 
data are 0(e), then the breaking times are respectively 0(e~^), 0(e~^) and 
0(6"^ ), for n = 1,2 and 3. We remark that, from the knowledge of such 
a breaking longtime regime of the solution, one can reconstruct exactly the 
initial data, an important issue in many physical contexts. 

The existence of a critical dimensionality above which small and localized 
data do not break has a clear origin, since, in the model, two terms act in 
opposite way: the nonlinearity is responsible for the steepening of the profile, 
while the n — 1 diffraction channels, represented by the transversal Laplacian, 
have an opposite effect; for n = 1,2,3 the nonlinearity prevails and wave 
breaking takes place (but at longer and longer time scales, as n increases), 
while, for n > 4, the number of transversal diffraction channels is enough to 
prevent such phenomenon, for small data, in the longtime regime. 

The paper is organized as follows. In §2 we derive the exact solutions of 
the model and in §3 we use them to build a uniform approximation of the 
solution of the Cauchy problem of ([T]), for small and localized initial data, 
establishing that wave breaking takes place only if n = 1, 2, 3. At last, in §4, 
we discuss in great detail the analytic aspects of such a wave breaking. 

2 Exact solutions of the dKPn equation 

The universal properties of the dKPn equation, discussed in the Introduction, 
suggest its invariance under motions on the associated paraboloid. Indeed, 
it is easy to show that equation ([1]) admits the following Lie point symmetry 
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group of transformations 

n-l 

yj = yj-'^^j^^ j = 1, • • • - 1, 

where the (5j's are the arbitrary parameters of the group, leaving invariant 
the paraboloid 

^ n— 1 

^+i^E^' = ^- (12) 

1=1 

Correspondingly, equation ([1]) possesses the following exact solution 

u={ A n-r' { (13) 

where F is an arbitrary function of one argument; such solution is charac- 
terized by the differential constraint a = 0, where a is the corresponding 
"characteristic symmetry" 



a 



'n~\ \ n-l 

^iyi I - 2t ^ 5iUy^ (14) 
A=l J i=l 

of equation ([T]). Indeed, if one looks for solutions of ([T]) in the form 

n-l 

At 



u = v{u), ^ = ^ + YiT.yi^ (15) 

i=l 

one obtains the following equation for v{^,t): 

11 — 1 

vt H —V + vvj: = 0. (16) 

Its v/t term can be eliminated by the change of variables 

v{^,t) = t-^q{^,T), (17) 

where ^ 

r(t) = I ^^3, .^g. 

^ ^ \ Int, n = 3, ^ ^ 
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leading to the Riemann - Hopf equation: 



Qr + qq^ = 0, 



(19) 



whose general solution is implicitely given by 



g = F(e-gr), 



(20) 



where F is an arbitrary function of one argument. Going back to the original 
variables via (HID, ([IB]), the solution becomes ([T^. 

We remark that, if F is a regular and localized function of its argument, 
the solution f|T3|) describes a wave concentrated on the wave front, given by 
the paraboloid f|T2l) . and constant on it, breaking, simultaneously, on the 
whole paraboloid. We also remark that, for n = 2, the solution f[T51) has been 
first derived using a recently developed method to construct exact solutions 
of the nonlinear RH inverse problem associated with the integrable dKP2 
equation [T9] . 

3 The Cauchy problem and wave breaking of 
small and localized initial data 

Since the paraboloid f|T2l) plays an important role in the asymptotics of the 
dKPn equation (see the Introduction), the exact solution f[T^ is physically 
relevant and can be used to build a uniform approximation of the solution of 
the Cauchy problem for the dKPn equation, under the hypothesis of small 
and localized initial data. 

The basic idea (already used in [20], in a different context) is that, if the 
initial condition is small: 



the solution of the Cauchy problem for dKPn is well approximated by the 
corresponding solution for the linearized dKPn untill one enters the nonlin- 
ear regime, in which the Riemann - Hopf equation f|T9l) becomes relevant, 
eventually causing wave breaking. Since the breaking time of 0(e) initial 
data evolving according to f|T9l) is r(t) = 0(e^^), the nonlinear regime for 
dKPn is characterized by the condition t = 0(r^^(e^^)), where is the 



u{x, y, 0) = euo{x, y), < e << 1, 



(21) 
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inverse of (fTSj) : so that: 

I e II n = 3, 

and a proper matching has to be made between the solution of the hnearized 
dKPn, vahd for t << 0(r~^(e~^)), and the exact solution of the previous 
section, valid in the nonlinear regime t = 0(r~^(e~^)). 

3.1 The linear regime 

Since the initial condition fl^T]) is small, the solution of dKPn is well approx- 
imated, for finite times, by the solution of the linearized dKPn equation: 

u{x,y,t) [ ^o(A:i,fc±)e^^'=^^+'^"-^^-^*)rfA;irffc^ (23) 

where Uo{ki,k±) is the Fourier transform of Uo{x,y) 

uo{kiM)= f Uo{x,y)e-'^^'''^^^-y^dxdy. (24) 

Such approximation is also valid in the longtime interval: 

1 « t « O (r-i (e-i)) (25) 

(far away from the nonlinear regime), in which the solution of dKPn is de- 
scribed by the standard stationary phase approximation of the multiple in- 
tegral ([23]): 

w(x, f, t) ~ r VeG (^x + l^ E yl l) , (26) 

where 

G(e,r/) :=2-"7r-^/^dA|A|'^Uo(A,Ar/)e^^«-^?("-i)^is"\ (27) 

valid in the space-time region ( 125|) and 

{x-O/t, y^|t = 0{l), z = l,...,n, (28) 

on the paraboloid f|T2l) . Outside the paraboloid, the solution decays faster. 
This formula says that the localized initial condition fj2T|) has evolved con- 
centrating, asymptotically, on the paraboloid ffT^ . 
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3.2 The nonlinear regime 

The approximate solution of dKPn in the nonhnear regime t = 0(r^^(e^^)), 
obtained matching equations fl26|) and f|T3l) . reads as follows: 



Since the term (ut) inside the first argument of function G, responsible 
for the wave breaking, is 0{t^~) for n 7^ 3, and the analogous term utlnt is 
O(lnt) for n = 3, then these terms are large for n = 1, 2, 3 and infinitesimal 
for n > A. It follows that wave breaking takes place only for n = 1,2,3; 
for n > 4 the solution (12^ coincides with the linearized solution fl2^ . and 
no breaking takes place. We remark that, in the integrable case n = 2, one 
recovers the results obtained in [8] using the 1ST for vector fields. We also 
remark that, if the breaking regime ( 129|) is known (measured), i.e., if function 
G is known, the initial condition euo{x,y) is uniquely reconstructed simply 
inverting (1271) . 

It is possible to show that the error made approximating the solution of 
dKPn by (EH]) is given by u = u'^^x, y, t){l + 0{t-^)) for n = 2, 3. 

Summarizing, the asymptotic solution fl29l) illustrates the following break- 
ing picture for the dKPn equation ([1]), corresponding to localized and 0(e) 
initial data. If n = 1 (the Riemann - Hopf case), waves break in the longtime 
regime t = 0(e~^); if ra = 2, waves break in the longtime regime t = 0(e~^), 
much later than in the 1 + 1 dimensional case; also if n = 3 small waves break, 
but at an exponentially large time scale: t = 0{e'^); at last, if n > 4, small 
and localized initial data do not break in the longtime regime. This result 
has a clear physical meaning: increasing the dimensionality of the transversal 
space, the number of diffraction channels of the wave increases, untill such 
diffraction, acting for a long time, is strong enough to prevent the gradient 
catastrophe of the small n dimensional wave. It is a remarkable coincidence 
that small initial data break, in the longtime regime, only in 1 + 1, 2 + 1 and 
3 + 1 dimensions; i.e., only in physical space! 

We end this section remarking that, if the initial condition is the Gaussian 
uo{x,y) = dn exp{— ^ ), where dn is a constant, then the above asymp- 
totic solution can be written in terms of elementary or special functions, 




(29) 
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depending on n: 



v^(l+|r7P)- 



COS^^i [ — ) 1^1 (^ — ,21 



y2 
4 



sm 



7r(ra-l) p (n+3\ y F f Sd^ 3- 



(30) 



where F is the Euler Gamma function and iFi is the Kummer confluent 
hypergeometric function [21]. If, in particular, n = 2 and n = 3, with 
(i2 = V2tt and = 2, one obtains, respectively. 



G(e,r/) = (l + r/2)-! 
(see Figure 1) and 



r (!) li^. 



3 1. _yf 

1 I 4' 2' 4 



+ 



(31) 



(l+r;2+r,2)^ 



4(l + r,^+r,|) _ 



(32) 




Figure 1. The plotting of function ( 1311) . 
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4 Geometric and analytic aspects of the wave 
breaking 

In this section we show how to derive, from the asymptotics fl29l) . the analytic 
features of the wave breaking for n = 2,3, in terms of the initial data, 
represented by function G defined in ( HT^ . We first rewrite equations fl2^ in 
the more convenient form 



w ~ eG{C,V), 



(33) 



where 



for n = 2, and 



w = Viu, r = 2Vi, ^ = x + ^^y', = ^ (34) 



w = tu, T = \nt, ^ = X + ^{yf + y^), V = (35) 

for n = 3, describing the evolution of an n- dimensional wave according to the 
1 + 1 dimensional Riemann - Hopf equation Wr + ww^ = 0. In the following, 
we mainly concentrate on the case n = 3; the case n = 2, that can be easily 
recovered setting to zero all the partial derivatives of G with respect to 7]2 in 
the formulae of this section, has been discussed in detail in [8]. 

One solves the second of equations f p3|) with respect to the parameter 
C, obtaining Ci^^V^''')^ and replaces it into the first, to obtain the solution 
w ~ eG{({^, f], T),f]). The inversion of the second of equations fl55]) is possible 
iff its ^-derivative is different from zero. Therefore the singularity manifold 
(SM) of the solution is the n - dimensional manifold characterized by the 
equation 

5(C,r/,r) = l + eGc(C,r/)r = => r = - ] (36) 

Since 

Vk*- = , (37) 

the slope of the localized wave becomes infinity (the so-called gradient catas- 
trophe) on the SM, and the n-dimensional wave "breaks" . 
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Then the first breaking time Tb and the corresponding characteristic pa- 
rameters = {(b, Vb) are defined by the global minimum of a function of n 
variables: 

n = ^ = global min (- ^ A > 0, (38) 

and it is characterized, together with the condition Gi^{(b) < 0, by the condi- 
tion that the symmetric quadratic form < z, Hz > be positive V2; G M", where 
"H is the Hessian matrix of function Gi^{(,f]), evaluated at C = Cfe = {(b,Vb)- 
The corresponding point at which the first wave breaking takes place is, 
from dMD, 6 = (6, ffb) G K", where: 

^b = Cb + eGiCbW (39) 

Now we evaluate equations (l33l) and (136!) near breaking, in the regime: 

e = 6 + r, V = Vb + V, r = Tb + T', C = Cb + C, (40) 



where ^' , rf , t' , (' are small. Using ( 138|) - (!39l) . the second of equations (133 
becomes, at the leading order, the following cubic equation in 



+ a{ff)C + K^, f)C - lX{i\ ff, f ) = 0, (41) 

where 

Kv', ^) = [2^^^ + G^r,ir,ir/'i^ + 2G^r^,^^r][r]!2 + Gc_r^^^^rl^\ , 
^(C', ^, r) = - eG(a, + rr)r' - e [G(a, + rf ) - G] Tfe ~ 



and 



r' r - Tb 



corresponding to the maximal balance 



(43) 



\aWMri'2\ = 0{n''), 1X1=0(1?^). (44) 
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In (1421) and in the rest of this section, all partial derivatives of G whose 
arguments are not indicated are meant to be evaluated at (h = {Cb, Vb)- 

The three roots of the cubic are given by the well-known Cardano-Tartaglia 
formula: 



C'i(r,rf,r) = -| + (A+)l + (A_)i, 



where 



A± = R±^/A 



and the discriminant A reads 



(45) 

(46) 
(47) 



with 



36- 



2|Gcl~ , 1 
7^ ' T 



G. 



ccc 



CCC L 



R{C, ff, f) = iX(e', rf , f ) + ft + IQ{ff, f) 
At the same order, function S in fl36l) becomes 



(48) 



(49) 



Known (' as function of (^', ^, r) from the cubic (14T!) . the solution w and its 
gradient are then approximated, near breaking, by the formulae: 



w{^,f],T)r^eG{Cb + C',Vb + v'), 



w ~ e 



(50) 



4.0.1 Before breaking 

If r < Tf, (f < 0), the coefficient Q in f HHj) is strictly positive, due to the 
positivity of the Hessian quadratic form; then the discriminant A = + 
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is also strictly positive and only the root d is real. Correspondingly, the real 
solution w is single valued and described by Cardano's formula (see Figure 
2). In addition, function S in fH9|) is also strictly positive and V(g,,f)W is finite 

To have a more explicit solution, we first restrict the asymptotic region 
to a narrower volume, so that the cubic ( 1411) reduces to the linear equation 
b(' = 7X; then the solution exhibits a universal behaviour, coinciding with 
the following exact similarity solution of equation Wr + ww^ = 0: 



describing the hyperplane tangent to the wave , where 



(51) 



h^,^) (52) 



defines the breaking direction. In addition, in such a narrow volume: 

'^{i,rf)W ^ — z7. (53) 

We also look at the different balance: rj'j of the same order, and f < 0(|r7^|), 
suitable for taking the r f Tf, limit. In this case the cubic simplifies to 
C'^ ~ 7X and 



G{c, + ^/WW^W^,Vb + v') 



(54) 



4.0.2 At breaking 



As r t Tfe, the above tangent hyperplane (now tangent at the breaking point) 
has an infinite slope and equation G^^' + Grfj^ri[ + Gr,2V'2 = 0- 

At the breaking time t = one can give an explicit description of the 
vertical inflection. If l^'/Vil; W/v-il = the cubic (HTl) simplifies to 



('^ = 7X and equation ( 15^ becomes 

w^eG (a + Vb + rf) ^ V^,,^ ~ ¥^ (55) 

V -^b 

where Xb{^',ff) = X{^',ff,0) is defined in ({421). We remark that Xb ~ 
r + i^Vi + §f^2), if C = c^iV'i + c^2V2, ^ -G,,/G^, J = 1,2, while 
~ ^{GviviV'i^ + '^GviV2VW2 + <^»72'?2'72^) if = -Gr,JGt^, j = 1, 2. 
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Equation fl55l) implies that, if n = 2, all the derivatives of w blow up at 
T = Tfc, in the breaking point {^b^Vb), with the universal law X^"^^^, except 
the derivative along the transversal line Xf,{^',7]') = 0, represented by the 
vector field X = X^^d^ — drj, for which 

^w\(i;b,vb) = -eCr,. (56) 

If n = 3, the situation is similar: all derivatives of w blow up at r = Tf,, in the 
breaking point {Cb,Vb), with the universal law except the derivatives 

along the transversal surface Xb{^', ff) = 0, having as natural basis the vector 
fields Xj = Xbri^d^ ^ ^vj^ J = 2, for which: 

= -eCr,,, 3 = 1,2. (57) 



4.0.3 After breaking 

After breaking, the solution becomes three-valued in a compact region of 
the (^,57) - space (see Figures 2, 3, 4), and does not describe any physics; 
nevertheless a detailed study of the multivalued region is important, in view 
of a proper regularization of the model. 

If T > Tf, (f > 0), in the regime dB]), the SM equation S = 0: 



2\G(;\t - GcCcC^ + G(^r,ir,iV'l^ + C^r^a r?2 ^2 ^ + 



(58) 



describes an ellipsoidal paraboloid in the {(', ff, t) space, with minimum at 
the breaking point {^b^'^b)- 

Eliminating C,' from equations ( 158|) and (HI]), one obtains the SM equation 
in space-time coordinates, coinciding with the A = condition, where A, Q 
and R are defined in (147|) and (HHj) . 

For n = 2, the SM is a closed curve with two cusps in the (^, rf) - plane 
at the points 

±./J^^|P%^(t^,1)V? (59) 

(see Figure 3), corresponding to the conditions Q = R = 0,on which the three 
real solutions of the cubic coincide. On the remaining part of the closed curve 
A = 0, two of the three real branches coincide: wi ~ eG{(b + CijV)^ = 
W- ~ ^G{C,b + Outside the closed curve, A > and the real solution 
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w is single valued; inside the closed surface A < and the real solution w 
is three valued. We remark that the transversal and longitudinal widths of 

1 3 

such a closed curve are respectively 0{t^) and 0(r2). Therefore this curve 
develops, at r = Tfe, from the breaking point {C,b,Vb)y with an infinite speed 
in the tranversal direction, and with zero speed in the longitudinal direction, 
recovering the results obtained in [S]. 

For n = 3, the SM is a closed surface in the (CyV) - space made of two 
surfaces having the same boundary: the transversal closed curve Q = R = 
(see Figure 4), on which the three real solutions of the cubic coincide. 
The Q = condition defines an ellipse in the {rii,ri2) plane with semi- axes 
2|G'clG'ccc , where 



A 



r, 



± J '~ I ' , _ _ (60) 



As before, in the remaining part of the closed curve A = 0, two of the three 
real branches coincide. While the axes of the transversal closed curve Q = 
R = are of 0(f2), the thickness of the longitudinal region between the two 

3 

surfaces is of 0(f2). Therefore this closed surface develops, at r = Th, from 
the breaking point {C,b, fjb), with an infinite speed in the tranversal directions, 
and with zero speed in the longitudinal direction. Intersecting this closed 
surface with any plane containing the C, - axis, one obtains a closed curve 
with two cusps as in Figure 3; therefore the closed curve Q = R = is made 
of all these cusps. 

We end these considerations remarking that the similarity solution be- 
fore breaking, the vertical inflection at breaking, and the compact three- 
valued region after breaking make clear the universal character of the gradi- 
ent catastrophe of two- and three-dimensional waves evolving according to 
the Riemann-Hopf equation (see Figure 2). 

Since the transformations (!34|) .( l35|) are globally invertible, for t ^ 0: 

1 r7^r^ nr"^ 
n = 2: u= -^w{^,t],t), t = —, x = ^ —, V = —, (61) 

n = 3: u = jw{^,f],T), t = e^, x = - e'^ {r]l + r]l) , y = 2e'^r], (62) 

all the above considerations are easily transfered to the dKPn case. In par- 
ticular, small and localized initial data evolving according to the dKPn equa- 
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tion dU) break, at tb = in the point {xb,yb) = {^b - TbVb/'^jVbTi /2) of 
the parabohc wave front x + y'^j^Atb) = if n = 2 [8], and at tb = e^'' 
in the point {xb,yb) = {C,b — e'^H^bP, 2e'^''f/f,) of the paraboloidal wave front 
^ + (1/1+ vD/i^tb) = ^b if n = 3. In addition, all the previous considerations 
concerning the universal character of such a wave breaking: the similarity 
solution before breaking, the vertical inflection at breaking, and the compact 
three-valued space region after breaking, are transfered in a straightforward 
way to the dKPn equation, for n = 2, 3. 

For instance, at t = and in the space region \x' + rjif^y'i + ^^261/2 1 = 
0{\yj\/ {2tb)) (the transformed of |^'| = 0(|?7^|)) equations (!55|) become 



where x' = x — Xb, if = y — yb and Xb{x',if) = Xb{i\ff). Again, if n = 2, 
all derivatives of u at the breaking point (x^, yb) blow up, except that along 
the transversal line Xb{x' ■, y') = 0, for which: 



Analogously, if n = 3, all derivatives of u at the breaking point {xb,yb) blow 
up, except those along the transversal surface Xbix',]/) = 0, represented by 

the basis vector fields Xj = Xby .dx — dy,, j = 1, 2, for which: 




(63) 




(64) 



e 



2 ' j 1,2. 



(65) 



2tb 



17 




18 



idM - idM 




Fig. 3. For n = 2 and G{(,ri) given by (!3T]) . two consecutive snapshots, 
at f = 0.1 and f = 0.3, describing the evolution of the three- valued region, 
from the breaking point {^b,Vb) ~ (3.26,0) (the center of the figures). This 
region is delimited by a closed curve with two cusps. 




Figure 4. For n = 3 and G{(, rji, rj2) given by ( [32l) . two consecutive snapshots, 
at f = 0.3 and f = 0.5, describing the evolution of the compact three- valued 
region from the breaking point (the center of the figures) . The closed surface 
delimiting the three-valued region is made of two surfaces with the same 
boundary, the closed curve Q = i? = 0. 
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Fig. 5. For n = 2, e = 0.2 and G{(, rj) given by fl3T|) . a detail of the parabolic 
wave front of dKP2 at t = t^, around the breaking point [x^, yb). 
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